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Abstract 

Invariant mass distributions of and combinations are investigated 

in order to study excited B mesons. The analysis is based on a data sample 
corresponding to 3.0 fb“^ of pp collision data, recorded by the LHCb detector at 
centre-of-mass energies of 7 and 8TeV. Precise measurements of the masses and 
widths of the Hi(5721)^’+ and B2{57A7)^’~^ states are reported. Clear enhancements, 
particularly prominent at high pion transverse momentum, are seen over background 
in the mass range 5850-6000 MeV in both B^tt~ and B^tt~^ combinations. The 
structures are consistent with the presence of four excited B mesons, labelled 
Hj(5840)^’'’“ and Hj(5960)*^’'’', whose masses and widths are obtained under different 
hypotheses for their quantum numbers. 
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1 Introduction 


The properties of excited B mesons containing a light quark can be described in the context 
of heavy quark effective theory (HQET) [^. Since the mass of the b quark is much larger 
than the QCD scale, the Lagrangian can be expanded in powers of l/rrih, where the leading 
term defines the static limit (m^ —)■ cxd). In the heavy quark approximation, the B mesons 
are characterised by three quantum numbers: the orbital angular momentum L (S, P, D for 
L = 0,1, 2 respectively) of the light quark, its total angular momentum jg = |L ± ||, and 
the total angular momentum J = | jg ± 11 of the B meson. The spectroscopic notation has 
the form where S' = 0 or 1 is the sum of the quark spins and where the quantum 

number n describes the radial excitations of the state. The PDG notation [^| (which 


is nsed in this paper) has the form B^j\m) or B^P{nL)^ where m is the mass in nnits 
of MeVQthe * snperscript is given to those states with natnral spin-parity P = (—1)"^ 
(J^ = 0+, 1“, 2+,...), and the subscript J is omitted for pseudoscalar and vector states. A 
prime may be used to distinguish two states with the same quantum numbers. 

For L = 0, there are two possible (J; jq) combinations, both parity-odd, corresponding 
to the B meson gronnd state with = 0“ and to the excited B* state with = 1~. 
Higher excitations are collectively referred to as B** states and decay strongly to lighter 
B mesons and pions. For L = 1 there are four different possible (J; jq) combinations, all 
parity-even. Predictions for the masses of such states and higher excitations spread over 
a wide range of valnes, as shown in Fig. 0mii- As can be seen in Fig. the states 
come in donblets (two values of J for each jq), and within each donblet, one has natnral 
and one unnatural spin-parity quantum numbers. States with natural spin-parity (except 
for 0+) can decay to both Btt and B*7i final states. States with nnnatural spin-parity 
cannot decay to the psendoscalar-pseudoscalar Bti final state due to parity conservation, 
bnt may decay to B*7r (Table [^. Since the B* meson decays to B'j, the signature from 
a doublet of B** states is given by three peaks in the Bn mass spectrnm (nnless the 
doublet includes a 0^ state): one from the natural spin-parity state decay to Bn, and 
two from both states decaying to B*n with a missing photon. Due to the missing photon, 
the peaks from B*n decays are shifted down from the trne B** mass by the difference 
between the B* and B masses (this feature recently allowed a precise determination of the 
B* — B mass difference from the B^K~ spectrum [^). Depending on the widths of the 
states and the mass resolution, two or all three of these peaks may overlap and be hard to 
distingnish experimentally. The B^ and B[ states are predicted to be very broad [3, 10 


(*)/ 


since they decay via S-wave (the comparable states in the charm sector have widths of 
aronnd 300 MeV D- However, the Bi and i ?2 states decay only via D-wave and are 
predicted and observed to be much narrower. Higher states snch as the H(2S), 

H*(2S), i? 2 (lD) and i? 3 (lD) are predicted to have widths in the 100-200MeV range [lo] , 
consistent with the recent measurement of the properties of the D* 3 ( 1 D) state 


In contrast to the sitnation in the charm sector, there is relatively little experimental 
information concerning B meson spectroscopy. The Hi(5721)° and H2(5747)° states 


^Natural units where h = c = 1 are used. 
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Figure 1: Mass predictions of the excited B states [S- 10 . The boxes cover the range of predictions 
for the masses of each state, and the red dots indicate the measured values. The horizontal lines 
correspond to the Btt (red) and B*tt (blue) thresholds. 


Table 1: Allowed decay modes for the excited B states. 



Allowed decay mode 


Btt 

B*71 

0+ 

yes 

no 

0-,l+,2-,... 

no 

yes 

l-,2+,3-,... 

yes 

yes 


have been observed by the CDF 14 and DO experiments, and recently the CDF 
collaboration has presented results on the charged isospin partners, together with evidence 
for a higher mass resonance 16 . This result has prompted theoretical speculation about 
the origin of the new state |17f 21 . While in the D meson system amplitude analyses of 
excited states produced in B decays can be used to determine their spin and parity (see, 
for example, Refs. [^[^|^), in the B meson system it is very difficult to assign with 
certainty quantum numbers to observed states. The labelling of the states follows the 
quark-model expectations for the quantum numbers, which have not been experimentally 
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verified. 

In this paper, the results of a study of B^ti~ and combinations are presented. 

The inclusion of charge-conjugate processes is implied throughout. The analysis is based 
on a data sample corresponding to 3.0 fb“^ of LHC pp collision data recorded with the 
LHCb detector at centre-of-mass energies of 7 and 8 TeV. 

The B mesons are reconstructed in the 

and D~n~^7r~^7T~ channels, with subsequent J/'il) —)■ —)• iC+vr” and iL+Tr^vr+Tr", 

D~ —)■ and K*^ —)■ K^tt~ decays. The B meson candidates are required to 

originate from a primary pp collision vertex (PV), and are combined with pions originating 
from the same PV (referred to as “companion pions”). Both “right-sign” (RS) and “wrong- 
sign” (WS) combinations are considered, where the latter are those with quark-content 
that precludes that the pair originates from the strong decay of an excited B meson {e.g. 
B~^7 t~^) and are used to model the combinatorial background. Excited B mesons are seen as 
peaks in the RS invariant mass distributions, and are htted with relativistic Breit-Wigner 
(RBW) functions. An additional very broad component, observed in the RS and not in 
the WS combinations, is referred to as “associated production” (AP) in this paper. The 
AP contribution may originate from very broad resonances or from correlated nonresonant 
production of B mesons and companion pions in the fragmentation chain. 

The remainder of the paper is organised as follows. A brief description of the LHCb 
detector is given in Sec. The selection requirements are described in Sec. the ht model 
is discussed in Sec. and the nominal fit results are given in Sec. with the evaluation 
of the systematic uncertainties in Sec. Interpretation of the results and a summary are 
given in Sec. 


2 Detector and dataset 


The LHCb detector |23[|24| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r/ < 5, designed for the study of particles containing 6 or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 25 


surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes 26 placed downstream of the magnet. The tracking 


system provides a momentum measurement with relative uncertainty that varies from 
0.5% at low momentum to 1.0% at 200 GeV, and an impact parameter measurement 
with resolution of 20 pm for tracks with large momentum transverse to the beamline 
(Pt)- Different types of charged hadrons are distinguished using information from two 
ring-imaging Cherenkov detectors j^. Photon, electron and hadron candidates are 
identihed by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identihed by a 


system composed of alternating layers of iron and multiwire proportional chambers 28 


The trigger 29 consists of a hardware stage, based on information from the calorimeter 


and muon systems, followed by a software stage, which uses information from the vertex 
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detector and tracking system. 

In the simnlation, pp collisions are generated nsing Pythia 
confignration 31 . Decays of hadronic particles are described by EvtGen 
hnal-state radiation is generated using Photos I^. The interaction of the generated 


with a specihc LHCb 
in which 
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particles with the detector, and its response, are implemented using the Geant4 toolkit 
as described in Ref. 
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3 Event selection 


The B meson candidates in each decay mode are reconstructed using a set of loose selection 
requirements to suppress the majority of the combinatorial backgrounds. The selection 
criteria are similar to those used in previous analyses of the same channels 36 39 . The 


—)■ J/'ipK^ and B^ —)■ J/ipK*^ selections require a B candidate with px > 3GeV 
and a decay time of at least 0.3 ps. For the other decay modes, the selection explicitly 
requires that the software trigger decision is based only on tracks from which the B meson 
candidate is formed. No requirement is imposed on how the event was selected at the 
hardware trigger stage. Additional loose selection requirements are placed on variables 
related to the B meson production and decay, such as transverse momentum and quality 
of the track hts for the decay products, detachment of the B candidate from the PV, and 
whether the momentu m of the B candidate points back to the PV. Because B^ mesons 
oscillate, the distinction between RS and WS combinations is clearest at short B^ decay 
times, and hence only B^ candidates with decay time below 2 ps are used in the analysis. 

The mass distributions for the B^ and B^ candidates are shown in Fig. Only B meson 
candidates falling within 25 MeV of the nominal B mass for the decay modes containing 
J/'ip mesons, or 50 MeV for the other modes, are selected for further analysis. Samples of 
about 1.2 million and 2.5 million B^ candidates are obtained, with purity depending on 
decay mode and always larger than 80%. Each candidate is combined with any track that 
originates from the same PV and that is identihed as a pion. The particle identihcation 
requirements on the companion pion are chosen to reduce potential backgrounds from 
misidentihed particles to a level where they can be neglected in the analysis. Over the 
momentum range relevant for this analysis, the pion identihcation requirements are 81% 
efficient at identifying pions, while they have 3.1% and 2.6% probabilities respectively to 
misidentify a kaon or a proton as a pion. Since the production of mesons is likely 
to be suppressed relative to the production of B** states, as has been observed for the 

these requirements are expected to reduce background from the 
—>■ B*^K~ and 5*9(5840)° — B*^K~ or B^K~, where the kaon is 


ground states |41| 42 
decays 5si(5830)° —r xj jv cui<a j ->^2 
misidentihed as a pion, to a negligible level. 

Further selection requirements are placed on the B** candidate. The invariant mass 
and x^/ndf (ndf is the number of degrees of freedom) of the B** candidate vertex ht are 
calculated constraining the B candidates and companion pion to originate from the PV, 
and also constraining the known B meson mass, and the masses of intermediate J/V’, 
5° and D~ mesons in the B decay. The y^/ndf of the B** candidate vertex ht is then 
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Figure 2: Mass distributions of the and candidates reconstructed through (a) B'^ 
(7r~^, 7T~^Tr~^7T~), (b) B~^ —)• (c) B^ —)• D“(7r'’', 7r'’'7r'’'7r“), and (d) B^ —)• J/iIjK*^ 

decays. The J/if) , and D~ masses are constrained to their world average values [^. Results 
of fits are superimposed for illustration. The signal (dot-dashed red line) is modelled with a 
double Crystal Ball distribution, while the background (dashed black line) is modelled with 
a second-order polynomial. The total fit is shown as a solid blue line. 


required to be below 3.5. In order to reduce combinatorial backgrounds, the PV associated 
with the B** candidate is required to have fewer than 75 charged particles associated with 
it. The angle 9 is required to satisfy cos6' > —0.5, where 9 is the angle between the pion in 
the Btt rest frame and the opposite direction of the boost vector from the Bn rest frame 
to the laboratory frame. 

Finally, the companion pion is required to have more than (0.5) 5 GeV of (transverse) 
momentum, while the B candidate is required to have pt > 10 GeV for candidates where 
the companion pion has pt > 2 GeV. In any selected event, the B candidate can potentially 
be combined with several different pions to create B** candidates. The average number of 
candidates per selected event is 1.4 and all of them are used for the subsequent analysis. 

4 Fit model 

The distributions of the mass difference, Q = m{Bn) — m{B) — m{n), following these 
selection requirements are shown in Fig. [^for both RS and WS B** candidates, where ms 
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Figure 3: Distributions of the Q values of the B** candidates after the selection for the (top) B~^ 
and (bottom) B^ candidates. The white histograms represent the RS combinations, while the 
overlaid shaded red histograms represent the WS combinations. The right hand plots are made 
after applying an additional requirement oi > 2 GeV on the companion pion. 


and are the known masses of the B meson and the pion |^. All B decay modes are 
combined in Fig. [^and in the snbseqnent analysis. Two narrow peaks are seen in both 
B^'k~ and B^tt^ mass difference distribntions, corresponding to the i?i(5721)‘^’’'‘ —)■ B*tt 
signal overlapping with the i?2(5747)°’+ — B*t: decay, and the i?2(5747)°’+ —)■ Bn decay. 
In addition, an excess of RS over WS combinations aronnd Q ~ 500 MeV is particularly 
prominent after requiring the companion pion to have pt > 2 GeV. This peak could result 
from a combination of two heavier B** resonances, consistent with the expectation that 
B** states come in doublets, as described in Sec. the structure is further analysed as 
described below. Furthermore, a comparison with the WS distributions shows a very broad 
excess of RS combinations lying under the resonances, corresponding to AP as discussed 
in Sec. [H 

The Q-value distributions of B~^n~ and B^n^ candidates are fitted independently to 
determine the masses and widths of the various resonant signals. In order to increase 
sensitivity to the parameters of the high mass states, the fits are performed in three bins 
of companion pion 0.5 < px < 1 GeV, 1 < px < 2 GeV and px > 2 GeV. The hts 
minimise the total of the Q-value distributions (in bins of width 1 MeV) simultaneously 
for the three companion pion px bins. 

The combinatorial background shape is obtained from WS combinations. It has been 
checked that the WS background consists of purely combinatorial background by studying 
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Btt combinations in which a B meson from one event is combined with a companion pion 
from another event; consistent shapes are found. The WS Q-value distributions are htted 
with piecewise-dehned, smooth polynomial (“spline”) functions. The shape is fixed in the 
subsequent £t to the RS distribution, but the yield is allowed to vary. 

Resonances are modelled with RBW lineshapes [^, given by 


Rrbw(^) 


r(m) 


( 1 ) 


where m is the Bn invariant mass (which is trivially related to the Q valne), mg is the 
mass value for the resonanc^and r(m) is the mass dependent width 


r(m) = T 


mo 
0 — 
m 


( 9(^) 

\q{mo) 


21+1 



( 2 ) 


In the latter eqnation Tq is the natnral width, q{m) is the B or n momentnm in the 
rest frame of the resonance and I is the orbital angnlar momentnm between the B and 
n mesons. The Blatt-Weisskopf form factors J^i 44,45 acconnt for the fact that the 


maximnm angnlar momentnm is limited by the phase-space in the decay. Dehning the 
dimensionless quantity z{m) = q^{m)R^, where R is the effective radius, Ri is dehned as 


J^2 


1 z{mo) 
1 -I- z{m) 


{z{mo) - -F 92:(mo) 
{z{m) — 3)^ -I- 9z(m) 


(3) 


Depending on the ht model, the B** resonances are described by hve or six RBW 
shapes: 

• one for the Ri(5721)°’’'' —)■ B*n feed-down into the left narrow peak with width, 
yield, and mean free to vary in the fits; 

• one for the B2{^1^1)^'^ —?■ Bn signal (the right narrow peak) with width, yield, and 
mean free to vary in the hts; 

• one for the B2{5747)^'~^ —)■ B*n feed-down into the left narrow peak with width 
hxed to be the same as that of the B2{7i747)^'^ —?■ Bn signal, mean shifted from the 
B2{7)747)'^'^ —?■ Bn peak by the known B* — B mass difference, 45.0 ± 0.4 MeV [^, 
and relative yield in px bins constrained as described later; 

• two (or three) for the higher mass components, with widths, means, and yields 
free to vary in the hts (except in the three RBW case, where two of the means are 
constrained by the B* — B mass difference). 

^The mass difference mo — m{B) — m{n) is referred to as the mean fjL hereafter. 
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The alternative descriptions for the higher mass resonances are motivated by the lack 
of knowledge of their qnantnm nnmbers. As described in Sec. a donblet of states is 
expected to give rise to three peaks. For example, for the (S(2S), i?*(2S)) donblet the 
higher (lower) mass of the pair has natural (unnatural) spin-parity. The description with 
three RBW shapes, two of which are constrained to have means offset by the B* — B mass 
difference, is therefore a physically motivated choice, obtained by applying quark-model 
expectations to the new states. However, there are two possibilities for this configuration, 
since it may be either the lower or the higher of the states that gives rise to two peaks. 
The alternative, with only two RBW shapes, is an empirical model, that corresponds to 
the minimal choice necessary to obtain a satisfactory description of the data. This is taken 
as the default and is referred to hereafter as the empirical model, but results of alternative 
fits with three RBW shapes are also presented. 

The RBW shapes have several parameters which need to be hxed in the hts, in 
particular the spin and effective radius input to the Blatt-Weisskopf form factors. The 
Ri(5721)°’+ and i?2(5'^47)°’+ resonances are assigned spin 1 and 2, respectively, and are 
both assumed to decay via D-wave (/ = 2), while the two higher mass resonances are 
assigned spin 0 (/ = 0) in the default ht. The effective radius is hxed to 4 GeV“^ 

The mass resolution is around 2 MeV which is negligible compared to the natural widths 
(> 20 MeV) of the resonances, and is therefore not modelled. The variation of the signal 
reconstruction efficiency with Q value is described with a hfth-order polynomial function 
with parameters determined from simulation. All signal parameters except the yields are 
shared between the different bins and B meson decay modes, though the efficiency 
function is determined independently for each bin. 

The AP component is caused by correlations between the B meson and the companion 
pion, and as such is not present in either the WS sample or in a sample obtained by mixing 
B mesons and pions from different events. As there is no suitable data control sample from 
which it can be constrained, it must be empirically modelled. The AP is modelled by a 
sixth-order polynomial shape determined from simulation with an additional broad spin-0 
RBW function to account for possible data-simulation differences. The latter component 
is introduced since the modelling of fragmentation effects in the simulation is expected to 
be imprecise. 

The relative yields of B2{5747)^’~^ —)■ B*7i and Bti in each px bin are fixed according 
to the relative efficiencies found in simulation, so that the relative branching fraction 
ratios B{B2{57A7)^’~^ —)■ B*7r)/B{B2{57A7)^’~^ —)■ Bn) are free parameters of the hts. The 
WS and AP yields are freely varied in the hts, independently in each px bin. The RBW 
parameters of the AP shape are also freely varied; the remaining parameters are hxed to 
the values obtained from simulation to avoid instabilities in the hts. The ht procedure is 
validated using large ensembles of pseudoexperiments. 
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5 Fit results 


The results of the empirical model hts to the B** candidates integrated over the three 
Pt bins are shown in Fig. The resnlts are also shown split by pt bin in Fig. where 
the plots have been zoomed into the range below 800 MeV in order to emphasise the 
resonant structnres. The resnlts for the parameters of interest are reported in Table 
Note that the reported mean values correspond to the peak positions, and do not include 
any correction for the B* — B mass difference, but when a state is assumed to have natural 
spin-parity, and therefore gives two peaks, the mass value reported is that of the higher 
peak. The results are consistent for the charged and nentral states, as expected since 
the nncertainties are larger than isospin splitting effects. The results for the higher mass 
states depend on whether they are assumed to have natural or unnatural spin-parity, and 
the results with the alternative hypotheses are presented in Table For the purpose of 
labelling, and withont prejudice on their qnantum numbers, the lower of these states is 
referred to subsequently as the i?j(5840)°’’'' and the other as the i?j(5960)°’’'' state. 


Table 2: Results of the fits when two RBW functions are used for the Rj(5840)‘^’’'' and B j(5960)‘’’’'' 
states (empirical model). The mean p of each peak is given together with the width T and the 
yield Nstate- The parameters related to the AP and WS components are suppressed for brevity. 
All uncertainties are statistical only. Units of MeV for p and T are implied. 


Fit parameter 

B+n- 

S°7r+ 

5i(5721)0’+ /r 

263.9 ±0.7 

260.9 ±1.8 

Si (5721)°’+ F 

30.1 ± 1.5 

29.1 ±3.6 

S*(5747)°’+ /i 

320.6 ±0.4 

318.1 ±0.7 

S*(5747)°’+ F 

24.5 ± 1.0 

23.6 ±2.0 

( 5721 ) 0 .+ low Pt 

14200 ± 1400 

3140 ± 750 

(5721)0.+ mid pt 

16200 ± 1500 

4020 ± 890 

(5721)0.+ high pt 

4830 ± 470 

940 ± 260 

^55(5747)0.+ low Pt 

7450 ± 420 

1310 ± 180 

^52(5747)0.+ mid Pt 

7600 ± 340 

2070 ± 180 

^52(5747)0.+ high Pt 

1690 ± 130 

640 ± 80 

S(S*(5747)°’+ ^ S*7r)/S(S*(5747)°’+ ^ Bn) 

0.71 ±0.14 

1.0 ±0.5 

Sj(5840)°’+ /i 

444 ± 5 

431 ± 13 

Sj(5840)0’+ F 

127 ± 17 

224± 24 

Sj(5960)°’+ /i 

550.4 ±2.9 

545.8 ±4.1 

Sj(5960)°’+ F 

82 ± 8 

63 ± 15 

(5840)0.+ low Pt 

3200 ± 1300 

1630 ± 970 

(5840)0.+ mid Pt 

5600 ± 1000 

3230 ± 720 

(5840)0.+ high Pt 

3090 ± 550 

2280 ± 450 

(5960)0.+ low Pt 

3270 ± 660 

610 ± 240 

^Sj(5960)o.+ mid Pt 

4590 ± 610 

910 ± 250 

(5960)0.+ high Pt 

2400 ± 320 

500 ± 140 
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Figure 4: Result of the fits to the Q-value distributions for (top) and (bottom) 

candidates. The components are labelled in the legend. The normalised residuals (pulls) of the 
difference between the ht results and the data points, divided by their uncertainties, are shown 
underneath each plot. 
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Figure 5: Result of the fit to (left) B^tt~ and (right) B^Tr~^ candidates, split into (top to bottom) 
low, medium and high px bins, with ranges as labelled on the plots. The components are labelled 
in the legends. The fit pulls are shown underneath each plot. 
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Table 3: Results of the fits when the natural spin-parity hypothesis is assigned to (top quadruplet) 
the i3j(5840)'^’''' state or (bottom quadruplet) the Rj(5960)'^’''' state, so that three RBW shapes 
are used to model the broad resonances in the fit. The mean /x of each peak is given together 
with the width T. All uncertainties are statistical only. Units of MeV for /x and T are implied. 


Fit parameter 

B+n- 

B^n+ 

Sj(5840)‘^’+ /i 

471 ± 22 

455 ± 26 

Sj(5840)°’+ F 

107 ±20 

215 ±27 

5j(5960)°’+ /i 

552 ± 4 

547 ± 4 

5j(5960)°’+ F 

82± 10 

61 ± 15 

5j(5840)°’+ /i 

444 ± 5 

425 ± 15 

5j(5840)°’+ F 

119 ± 17 

229 ± 27 

5j(5960)°’+ /i 

575 ± 6 

547 ± 5 

5j(5960)°’+ F 

56 ± 7 

61 ± 14 


The covariance matrix of the empirical model ht is given in Appendix For brevity, 
the results for the signal yields and for the background parameters are not reported. The 
magnitudes of the correlations between the signal observables and background shapes are 
smaller than 30%. All hts have acceptable minimum values. 

6 Systematic uncertainties 

Systematic uncertainties are evaluated in a data-driven manner, by varying fit parameters 
or conhgurations from their default values and taking the difference in the ht results as a 
systematic uncertainty. Summaries of the systematic uncertainties are given in Tables 
and for B^tt~ and resonances.The total systematic uncertainties on each individual 

parameter are obtained by combining all sources in quadrature. The covariance matrix 
of the systematic uncertainties, given in Appendix is computed by considering the 
correlated effects on the ht parameters of the systematic uncertainties. 

The modelling of the background shapes may depend on the choice of ht range. The 
upper and lower edges of the range are varied independently by 20% to assign systematic 
uncertainties. Similarly, any dependence of the results on the choice of bin width is 
evaluated by repeating the hts with 2 (instead of 1) MeV binning. Additional uncertainties 
due to the background modelling are assigned by varying the spline function used to 
describe the WS distribution and by varying parameters of the AP polynomial function. 

The relative efficiencies of the i?2(5747) decays to B*tt and Bn in each of the three 
Pt bins are hxed in the nominal ht. These are varied independently to assign systematic 
uncertainties. The uncertainties in the dependence of the efficiency on Q value are 
propagated to the results by repeating the ht after varying, within their errors, the 
parameters of the polynomial function used to describe the variation. Uncertainties are 
assigned for possible differences between data and simulation in the efficiency function 
by reweighting the simulation to match the B momentum distributions observed in data. 
Uncertainties are also assigned to take in account the effect of changing the pt > 3 GeV 


12 






Table 4: Systematic uncertainties on the results of the fit to the candidates. Units of 

MeV for ^ and T are implied. 


Source 

Pi (5721)0 

P^(5747)° 


Pj(5840)° 

Pj(5960)° 


h 

r 

BF ratio 

h 

r 

h 

r 

h 

r 

Total statistical 

0.72 

1.52 

0.14 

0.37 

1.01 

4.95 

16.70 

2.88 

7.71 

Fit range (high) 

0.33 

1.30 

0.06 

0.08 

0.37 

2.20 

2.90 

0.52 

0.26 

Fit range (low) 

0.04 

0.11 

0.01 

0.02 

0.39 

0.04 

8.22 

0.69 

2.83 

2 MeV bins 

0.02 

0.14 

0.00 

0.04 

0.07 

1.09 

0.50 

0.08 

1.00 

Spline knots 

0.11 

0.01 

0.02 

0.02 

0.26 

1.75 

0.04 

0.45 

1.44 

Float AP 

0.03 

0.00 

0.00 

0.03 

0.30 

1.58 

10.16 

0.73 

4.23 

i?2(5747)° rel. eff., low px 

0.56 

0.91 

0.15 

0.08 

0.16 

0.07 

0.23 

0.00 

0.18 

i?2(5747)° rel. eff., mid px 

0.64 

1.01 

0.05 

0.09 

0.18 

0.08 

0.26 

0.00 

0.16 

i?2(5747)° rel. eff., high px 

0.20 

0.37 

0.03 

0.02 

0.07 

0.02 

0.00 

0.01 

0.09 

Eff. variation with Q value 

0.13 

0.33 

0.02 

0.04 

0.07 

0.45 

2.46 

0.19 

0.70 

Data-simulation reweighting 

0.07 

0.38 

0.02 

0.00 

0.16 

1.81 

2.03 

0.49 

0.12 

B Px 

0.02 

0.20 

0.01 

0.24 

0.72 

3.98 

3.67 

1.30 

4.29 

Px binning 

0.90 

2.45 

0.24 

0.06 

0.39 

1.49 

27.77 

4.20 

1.47 

Fit bias 

0.06 

0.17 

0.01 

0.00 

0.16 

0.45 

5.34 

0.40 

2.24 

Spin 

0.02 

0.06 

0.01 

0.06 

0.46 

1.95 

3.32 

0.62 

3.74 

Effective radius 

0.33 

1.44 

0.02 

0.12 

0.76 

2.17 

9.68 

1.24 

3.81 

B* — B mass 

0.10 

0.11 

0.03 

0.02 

0.11 

0.04 

0.17 

0.00 

0.09 

5j(5840)° 

0.01 

0.04 

0.00 

0.01 

0.01 

— 

— 

1.67 

0.76 

5j(5960)° 

0.01 

0.20 

0.00 

0.00 

0.16 

0.18 

8.00 

— 

— 

Extra state 

0.00 

0.26 

0.00 

0.04 

0.34 

1.67 

0.99 

0.12 

2.08 

Total systematic 

1.36 

3.49 

0.30 

0.33 

1.48 

6.68 

34.24 

5.10 

9.41 


cut on the B candidate to pt > 4 GeV, and of varying the boundaries of the three bins of 
the companion pion px- 

Possible biases in the hts are investigated with ensembles of pseudoexperiments. No 
signihcant bias is found for most of the parameters, but shifts in the means and widths of 
the i?j(5840)° and i?j(5960)° states of up to 30% of the statistical uncertainty are found 
and corrected for. Systematic uncertainties corresponding to the size of the bias seen in 
the ensembles are assigned to all parameters. 

Further systematic uncertainties are evaluated for the hxed £t parameters. The spins 
of the higher mass states are changed from zero to two, the Blatt-Weisskopf effective 
radius is varied from its nominal value of 4 GeV“^ to 2 and 6 GeV“^, and the B* — B mass 
difference is varied within its uncertainty [^. The effects on the other parameters of the 
ht, when the 5j (5840)° and 5j (5960)° states are assumed to have natural spin-parity 
and hence contribute two peaks to the spectrum, are assigned as systematic uncertainties; 
the effects on the parameters of the i?j(5840)° and i?j(5960)° states themselves when 
changing this assumption are presented in Table Finally, the hts are repeated allowing 
for an additional state with a peak around Q ~ 800 MeV. The additional state is not 
statistically signihcant, but the changes in the htted parameters are assigned as systematic 
uncertainties. The systematic uncertainties due to the momentum scale calibration are 
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Table 5: Systematic uncertainties on the results of the fit to the candidates. Units of MeV 

for and T are implied. 


Source 

5i(5721)+ 

5^(5747)+ 


5j(5840)+ 

Bj(5960)+ 


h 

F 

BF ratio 

h 

F 

h 

F 

h 

F 

Total statistical 

1.81 

3.57 

0.51 

0.72 

1.99 

12.70 

23.90 

4.07 

14.50 

Fit range (high) 

0.35 

0.74 

0.10 

0.11 

0.25 

1.51 

12.85 

0.38 

0.46 

Fit range (low) 

0.64 

1.13 

0.13 

0.06 

0.13 

7.85 

39.71 

0.14 

1.44 

2 MeV bins 

0.16 

0.34 

0.05 

0.10 

0.49 

0.58 

3.84 

0.28 

0.52 

Spline knots 

0.30 

0.08 

0.07 

0.03 

0.22 

1.94 

2.64 

0.25 

0.25 

Float AP 

0.02 

0.31 

0.01 

0.02 

0.03 

2.91 

2.44 

0.19 

2.24 

i?2(5747)+ rel. eff, low px 

1.50 

2.14 

0.43 

0.12 

0.49 

0.15 

1.63 

0.02 

0.03 

i?2(5747)+ rel. eff, mid px 

1.55 

2.26 

0.53 

0.12 

0.51 

0.29 

2.03 

0.04 

0.15 

132(5747)+ rel. eff, high px 

0.49 

0.90 

0.11 

0.03 

0.12 

0.10 

0.84 

0.02 

0.07 

Eff. variation with Q value 

0.07 

0.27 

0.02 

0.03 

0.10 

1.65 

7.28 

0.16 

0.94 

Data-simulation reweighting 

0.04 

0.38 

0.03 

0.00 

0.02 

2.13 

7.49 

0.40 

1.75 

B Px 

0.45 

1.38 

0.17 

0.14 

0.54 

1.16 

7.79 

0.98 

4.65 

Px binning 

1.82 

1.03 

0.26 

0.15 

1.38 

0.54 

55.56 

0.94 

11.43 

Fit bias 

0.14 

0.39 

0.04 

0.01 

0.32 

1.14 

7.65 

0.57 

4.21 

Spin 

0.14 

0.33 

0.05 

0.15 

0.94 

4.18 

24.49 

1.67 

5.98 

Effective radius 

0.70 

1.48 

0.12 

0.19 

0.29 

2.82 

22.15 

0.39 

3.76 

B* — B mass 

0.21 

0.06 

0.07 

0.01 

0.15 

0.32 

0.48 

0.03 

0.07 

5,7(5840)+ 

0.00 

0.05 

0.00 

0.03 

0.15 

— 

— 

0.72 

1.64 

5,7(5960)+ 

0.02 

0.01 

0.01 

0.04 

0.26 

5.99 

4.86 

— 

— 

Extra state 

0.03 

0.41 

0.00 

0.00 

0.15 

6.28 

12.82 

0.43 

7.81 

Total systematic 

3.10 

4.28 

0.79 

0.40 

2.07 

13.70 

79.82 

2.52 

17.18 


found to be negligible. 

In addition, varions cross-checks are performed to ensure fit stability and reliability. 
The stability of the data fits is studied by splitting the sample by the year of data taking, 
magnet polarity, and the charge of the companion pion. The resulting independent samples 
are fitted using the same conhguration as the nominal fit, and the results within each split 
are fonnd to be consistent. 


7 Interpretation and conclusions 

The analysis of the invariant mass spectra of B^tt~ and combinations reconstructed 

with the LHCb detector reported in this paper provides measurements of the properties of 
a nnmber of B** resonant states. The interpretation of the results is now given in two 
parts: firstly for the narrow B** signals, and subsequently for the broad, higher mass B** 
signals. 

The narrow states are identihed with the previously observed i?i(5721)° and i?2(5747)° 
states, and their i7i(5721)’'' and B2{5747)~^ isospin connterparts. The peak positions in 
the Q-valne distributions reported in Sec. [^can be converted into absolnte masses using 
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the known B and tt meson masses and the B* — B mass difference , leading to 


^Si(5721)0 

= 5727.7 

^B|(5747)0 

= 5739.44 

’^Si(5721)-I- 

= 5725.1 

(5747)4- 

= 5737.20 

Tbi (5721)0 

30.1 

rB|(5747)0 

24.5 

Tbi (5721)4- 

29.1 

rB|(5747)4- 

23.6 


± 

0.7 

± 

1.4 

± 

± 

0.37 

± 

0.33 

± 

± 

1.8 

± 

3.1 

± 

± 

0.72 

± 

0.40 

± 

± 

1.5 

± 

3.5 


± 

1.0 

± 

1.5 


± 

3.6 

± 

4.3 


± 

2.0 

± 

2.1 



0.17 ± 

0.4 

MeV, 

0.17 


MeV, 

0.17 ± 

0.4 

MeV, 

0.17 


MeV, 


MeV, 

MeV, 

MeV, 

MeV. 


The listed uncertainties are, from left to right: the statistical uncertainty, the experimental 
systematic uncertainty, and, where applicable, the uncertainty on the B meson mass and 
the uncertainty on the B* — B mass difference. Note that i?i(5721)°’+ and i?2(5747)°’’'' 
notations are maintained here for consistency with the previous literature, even though 
the values of the masses no longer agree with these labels within uncertainty. The results 
reported above are the most precise determinations of these quantities to date. 

The relative branching fractions for the B2{5747)^’~^ decays are measured to be 


i3(S*(5747)° ^ B*+7r-) 

(5^(5747)0 ^ B+tt-) 

B{B*{5747)+ B*^7r+) 

13 {B*2{5747)+ ^ B^7r+) 


= 0.71 ± 0.14 ± 0.30, 

= 1.0 ± 0.5 ± 0.8 , 


where the uncertainties are statistical and systematic, respectively. The significances 
of the i?2(5747)°’+ —>■ B*n decays are evaluated using a likelihood ratio test. Values of 
6.5(7 and 1.8(7 are obtained for B*~^7r~ and B*^7r~^, respectively, when only the statistical 
uncertainty is considered. The inclusion of systematic uncertainties reduces the signihcance 
for the B*^'k~ case to 3.7(7. This result therefore corresponds to the first evidence for the 
i?2(5747)° —j- B*^tt~ decay. The relative branching fractions for the i?2(5747)°’’'' decays 


are in agreement with theoretical predictions 10,46-49 


Structures at higher mass are clearly observed in the Q-value distributions. To 
investigate the significance of the high mass states, large samples of pseudoexperiments 
are generated and fitted with different configurations. To cover the dominant systematic 
uncertainty on the yield of these states which arises due to lack of knowledge of the shape 
of the AP component, the pseudoexperiments are generated with the AP shape that 
minimises the significance. A first ensemble is generated without any high mass states 
included. Each pseudoexperiment in this ensemble is fitted twice, once with the same 
model as used for generation and once with an additional high mass resonance included. 
The distribution of the difference of values between the two fits is extrapolated to 
obtain the p-value corresponding to the probability to find a difference as large or larger 
than that obtained from the corresponding hts to data. This procedure gives signihcances 
of 9.6(7 for the B^n~ case and 4.8a for the B^tt^ case. 
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A second ensemble of pseudoexperiments is generated with a configuration that corre¬ 
sponds to the best fit to the data with a single high mass resonance. The pseudoexperiments 
in this ensemble are htted both with the model used for generation and with a second 
high mass resonance included. The significances of the second peaks, again obtained from 
the difference in values, are found to be 7.5a and 4.6 ct for the B^'k~ and cases, 

respectively. Since isospin symmetry is expected to hold for these states, this shows that 
under the hypothesis that the high mass structures are due to resonances, two new pairs 
of particles are observed. 

Masses and widths of the i?j(5840)°’’'' and i?j(5960)°’’'' states are obtained with different 
ht models, as discussed in Sec. and the corresponding results are shown in Table 
The properties of the i?j(5960)states are consistent with and more precise than those 
obtained by the CDF collaboration when assuming decay to Bn 16 . If the i?j(5840) 
and i?j(5960)°’’'' states are considered under the quark model hypothesis, their properties 
are consistent with those expected for the i?(2S) and B*{2S) radially excited states. 

In summary, the B^n~ and B^n^ invariant mass distributions obtained from LHC pp 
collision data recorded at centre-of-mass energies of 7 and 8 TeV, corresponding to an 
integrated luminosity of 3.0 fb~^, have been investigated in order to study excited B mesons. 
Precise measurements of the masses and widths of the i?i(5721)°’+ and i?2(5747)°’’'' states 
are reported. Evidence is found for the B2{57A7)^ —)■ B*^n~ decay. Clear enhancements 
over background are observed in the mass range 5850-6000 MeV in both B^n~ and B^n~^ 
combinations. Fits to the data, accounting for the apparent enhanced production of the 
high mass states in the high transverse momentum region, allow the parameters of these 
states, labelled Hj(5840)°’’'' and Hj(5960)°’’'', to be determined under different hypotheses 
for their quantum numbers. 
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Table 6: Parameters of the Rj(5840)®’“’' and Rj(5960)*^’'’' states obtained with different fit models. 
The empirical fit uses two, and the quark model fits three, RBW shapes to model the broad 
resonances. The listed uncertainties are, from left to right: the statistical uncertainty, the 
experimental systematic uncertainty, and, where applicable, the uncertainty on the B meson 
mass and the uncertainty on the B* — B mass difference. Note that any state not explicitly 
labelled as “natural” is considered to have unnatural spin-parity (and not to be O"''); the reported 
mass can be converted into the corresponding result under the O"*" spin-parity assumption by 
subtracting the B* — B mass difference. Units of MeV are implied. 


Empirical model 


(5840)0 

5862.9 

± 

5.0 

± 

6.7 

± 

0.2 



Tbj (5840)0 

127.4 

± 

16.7 

± 

34.2 





^Bj(5960)0 

5969.2 

± 

2.9 

± 

5.1 

± 

0.2 



rBj(5960)0 

82.3 

± 

7.7 

± 

9.4 





^Bj(5840)+ 

5850.3 

± 

12.7 

± 

13.7 

± 

0.2 



rBj(5840)+ 

224.4 

± 

23.9 

± 

79.8 





^^Bj(5960)+ 

5964.9 

± 

4.1 

± 

2.5 

± 

0.2 



rBj(5960)+ 

63.0 

± 

14.5 

± 

17.2 






Quark model, ilj(5840) 

0,-r 

natural 


(5840)0 

5889.7 

± 

22.2 

± 

6.7 

± 

0.2 



Tbj (5840)0 

107.0 

± 

19.6 

± 

34.2 





^^Bj(5960)0 

6015.9 

± 

3.7 

± 

5.1 

± 

0.2 

± 

0.4 

rBj(5960)0 

81.6 

± 

9.9 

± 

9.4 





^^Bj(5840)+ 

5874.5 

± 

25.7 

± 

13.7 

± 

0.2 



rBj(5840)+ 

214.6 

± 

26.7 

± 

79.8 





^Bj(5960)+ 

6010.6 

± 

4.0 

± 

2.5 

± 

0.2 

± 

0.4 

rBj(5960)+ 

61.4 

± 

14.5 

± 

17.2 






Quark model, ilj(5960) 

o,+ 

natural 


^Bj (5840)0 

5907.8 

± 

4.7 

± 

6.7 

± 

0.2 

± 

0.4 

Tbj (5840)0 

119.4 

± 

17.2 

± 

34.2 





^Bj(5960)0 

5993.6 

± 

6.4 

± 

5.1 

± 

0.2 



rBj(5960)0 

55.9 

± 

6.6 

± 

9.4 





^Bj(5840)+ 

5889.3 

± 

15.0 

± 

13.7 

± 

0.2 

± 

0.4 

rBj(5840)+ 

229.3 

± 

26.9 

± 

79.8 





^Bj(5960)+ 

5966.4 

± 

4.5 

± 

2.5 

± 

0.2 



rBj(5960)+ 

60.8 

± 

14.0 

± 

17.2 
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A Covariance matrices 


Tables and 1^ each show both statistical and systematic correlations between the main 
parameters of interest in the B^ti~ and fits, respectively. In each Table, the masses 

and widths of the two broad states are seen to be heavily correlated with each other 
becanse they overlap, while the parameters of the narrow states are correlated becanse of 
the overlap between the i?i(5721)°’+ state and the i?2(5747)°’+ feed-down. 


Table 7: Statistical (top) and systematic (bottom) covariance matrices of the nominal B~^7r~ fit, 
where /U and T stand for the mean and width respectively. The parameters related to the AP 
and WS shapes and the signal yields are suppressed for brevity. Units of MeV for /r and T are 
implied. 



5i(5721)° 

B*(5747)° 


Bj(5840)° 

Bj(5960)° 


h 

F 

BF ratio 

h 

F 

h 

F 

h r 

5 i ( 5721 )° /r 

0.5 








5 i ( 5721)0 r 

0.8 

2.3 







5* (5747)° BF ratio 

- 0.1 

- 0.1 

0.0 






52*(5747)° 

0.1 

0.1 

0.0 

0.1 





5^(5747)° r 

- 0.2 

-0.4 

0.0 

0.0 

1.0 




5j(5840)° n 

0.0 

-0.4 

0.0 

0.0 

0.0 

24.5 



5j(5840)° r 

- 0.1 

2.0 

0.0 

-0.4 

- 1.2 

23.1 

278.9 


5j(5960)° n 

0.0 

0.1 

0.0 

0.0 

0.0 

7.4 

21.2 

8.3 

5j(5960)° r 

0.1 

0.6 

0.0 

0.1 

0.9 

-21.4 

-41.2 

-10.2 59.4 

“ 5 ^( 572177 ^ 

1.9 








5 i ( 5721 )° r 

1.0 

12.2 







5^(5747)° BF ratio 

- 0.1 

- 0.1 

0.1 






52*(5747)° 

0.1 

0.1 

0.0 

0.1 





52*(5747)° F 

- 0.2 

-0.3 

0.0 

0.0 

2.2 




5j(5840)° n 

0.1 

0.1 

0.0 

0.0 

0.0 

44.6 



5j(5840)° F 

- 0.2 

-0.4 

0.0 

0.0 

0.1 

0.0 

1172 


5j(5960)° n 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

26.0 

5j(5960)° F 

0.2 

0.3 

0.0 

0.0 

0.0 

0.0 

- 0.1 

0.0 88.6 
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Table 8: Statistical (top) and systematic (bottom) covariance matrices of the nominal fit, 

where fx and T stand for the mean and width respectively. The parameters related to the AP 
and WS shapes and the signal yields are suppressed for brevity. Units of MeV for n and T are 
implied. 



5i(5721)+ 

. 62 * (5747)+ 


5j(5840)+ 

5j(5960)+ 



F 

BF ratio 


F 

fJ- 

F 

^x F 

5i(5721)+ 

3.3 








5i(5721)+ r 

5.0 

12.7 







i?2(5747)+ BF ratio 

-0.9 

-1.5 

0.3 






i?2*(5747)+ fi 

0.4 

0.4 

- 0.1 

0.5 





5* (5747)+ r 

- 0.8 

-1.9 

0.2 

0.1 

4.0 




5j(5840)+ fj, 

0.5 

-3.2 

- 0.1 

1.6 

8.8 

161.3 



5j(5840)+ r 

2.2 

9.4 

-0.7 

-0.9 

-7.6 

-42.5 

571.2 


5j(5960)+ 

0.1 

0.0 

0.0 

0.2 

1.0 

20.7 

- 7.8 

16.6 

5j(5960)+ r 

-0.3 

1.0 

0.0 

-0.4 

- 2.0 

-95.8 

-107.4 

-22.4 210.2 

5i(5721)+ 

9.6 








5i(5721)+ r 

3.7 

18.3 







i?2(5747)+ BF ratio 

- 0.8 

- 1.1 

0.6 






i?2*(5747)+ 

0.2 

0.3 

- 0.1 

0.2 





B;{574:7)+ F 

- 0.8 

- 1.2 

0.2 

- 0.1 

4.3 




5j(5840)+ 

- 0.2 

-0.3 

0.0 

0.0 

0.1 

187.7 



5j(5840)+ F 

3.0 

4.3 

-0.9 

0.2 

- 1.0 

-0.3 

6371 


5j(5960)+ fi 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

6.4 

5j(5960)+ F 

0.0 

0.0 

0.0 

0.0 

0.0 

- 0.1 

0.0 

0.0 295.2 
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